Logarithmically Improved Blow-up Criteria for a Phase Field 
Navier-Stokes Vesicle-Fluid Interaction Model* 

Jihong Zhao a t and Qiao Liu b 

a College of Science, Northwest A&F University, Yangling, 
Shaanxi 712100, PR China 

b Department of Mathematics, Hunan Normal University, Changsha, 
Hunan 410081, PR China 

Abstract 

In this paper, we study a hydrodynamical system modeling the deformation of vesicle 
membrane under external incompressible viscous flow fields. The system is in the Eulerian 
formulation and is governed by the coupling of the incompressible Navier-Stokes equations 
with a phase field equation. In the three dimensional case, we establish two logarithmically 
improved blow-up criteria for local smooth solutions of this system in terms of the vorticity 
field only in the homogeneous Besov spaces. 

Keywords: Phase field; Navier-Stokes equations; fluid vesicle interaction; regularity cri- 
terion; Besov space 

2010 AMS Subject Classification: 35B44, 35Q30, 76D09, 76T10 

1 Introduction 

Recently, there have been many numerical and theoretical studies on the configurations and de- 
formations of elastic vesicle membranes under external flow fields [3J |H |51 HI M El UHl I2H Hi] . 
The single component vesicle membranes are possibly the simplest models for the biological cells 
and molecules and have widely studied in biology, biophysics and bioengineering. Such vesicle 
membranes can be formed by certain amphiphilic molecules assembled in water to build bilayers, 
and exhibit a rich set of geometric structures in various mechanical, physical and biological envi- 
ronment [7] . In order to model and understand the formation and dynamics of vesicle membranes 
and the fluid structure interaction, one approach is to consider equations of elasticity for the vesicle 
membranes and the Navier-Stokes equations for the fluid. However, the model established in this 
approach is very difficult to study and numerically simulate due to the fact that the description 
for elasticity is in Lagrangian coordinate (Hooke's law) and for fluids is in Eulerian coordinate. To 
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overcome this difficulty, in [H[7], the authors established a phase field Navier-Stokes vesicle fluid 
interaction model for the vesicle shape dynamics in flow fields via the phase field approach. In this 
model, the vesicle membrane T is described by a phase function <fi, which is a labeling function 
defined on computational domain Q. The function <\> takes value +1 inside of the vesicle membrane 
and — 1 outside, with a thin transition layer of width characterized by a small (compared to the 
vesicle size) positive parameter e. Obviously, the sharp transition layer of the phase function gives 
a diffusive interface description of the vesicle membrane T, which is recovered by the zero level set 
{x : 4>{x) — 0}. The advantage of introducing such a phase function is to formulate the original 
Lagrangian description of the membrane evolution in the Eulerian coordinates. On the other hand, 
the viscous fluid is modeled by the incompressible Navier-Stokes equations with unit density and 
with an external force defined in terms of 4>. 

In this paper, we study the three dimensional phase field Navier-Stokes vesicle-fluid interaction 
model subjecting to the periodic boundary conditions (i.e., in torus T 3 ), which reads as follows: 

d t u + u ■ Vu + VP = fiAu + — 7y^-V</> inQx[0,T], (1.1) 



V • u = in Q x [0,T], (1.2) 

E{<j> 
5(f) 



d t cj> + u ■ V(f> = -J*^P- inQx[0,T] (1.3) 



with the initial conditions 

u(x,0) = uq(x) with V • Mo = 0, and <f>{x, 0) = <fio(x) for x G Q, (1-4) 

and the boundary conditions 

u(x + e i: t) = u(x,t), 4>(x + ei,t) = (f>(x,t) for x G dQ x [0, T], i = 1, 2, 3, (1.5) 

where the set of vectors {e\ = (l,0,0),e2 = (0, 1,0), e3 = (0,0,1)} denotes an orthonormal basis 
of M 3 and Q is the unit square in R 3 . Here u = (ui, 112, 1*3) G R 3 and P — P(x, t) el denote the 
unknown velocity vector field and unknown pressure of the fluid, respectively. G R is the phase 
function of the vesicle membrane T. E(<p) denotes the physical approximation/regularization of 
the Hclfrich elastic bending energy for the vesicle membrane which is given by (cf . [4j [6j HI |9] ) 

E(<t>) = EM + \Mi{A{cj>) - a) 2 + hf 2 (B(0) - (3) 2 (1.6) 

with 

EM = 4- f l/WI 2 ^ and ftt) = ~£&4>+-{4> 2 - 1)6 (1.7) 
2e Jn £ 

where e is a small (compared to the vesicle size) positive parameter that characterizes the thickness 

of transition layer of the phase function, M\ and M% are two penalty constants which arc introduced 

in order to enforce the volume 

A{4>) = [ 4>dx (1.8) 

and the surface area 
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B{$) = J (||V0| 2 + l(0 2 -l) 2 )dx (1.9) 
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of the vesicle conserved (in time), and a — A((f>o) and /3 — B(<j>o) are determined by the initial value 
of the phase function <pQ. The positive constants v, k, and 7 denote, respectively, the viscosity of 
the fluid, the bending modulus of the vesicle, and the mobility coefficient. 5E ^ is the so-called 
chemical potential that denotes the variational derivative of E(<j>) in the variable <j>. Note that, if 
we denote 

ff ( ( /») = -A/( ( /)) + l(30 2 -l)/(0), (1.10) 

e z 

then a direct calculation yields that the variation of the approximate elastic energy is given by (see 

Hi) 

5 -^- = kg{4>) + Mt(A((f>) ~a)+ M 2 (2?(0) - /?)/(</>) 

= fceA 2 - -A(0 3 - <f>) - -(30 2 - 1)A0 + 4(3^ 2 - 1)(0 2 - l)<f> 
+ M X {A{4>) - a) + M 2 (B(<t>) - p)f{4>). (1.11) 

The system (|1.1|> — (|1.3|l describes the dynamic evolution of vesicle membranes immersed in an 
incompressible, Newtonian fluid, using an energetic variational approach [H [7] (see El El EH 
[24] for numerical simulations and other studies). Equations (|1.1|) and ()1.2|) are the momentum 
conservation and the mass conservation equations of a viscous fluid with unit density and with an 
external force caused by the phase field <f>. Equation (|1.2p is the condition of incompressibility. 
Equation (( 1 . 3|) is a relaxed transport equation of 4> with advection by the velocity field u. The 
right-hand side of (|1.3j) is a regularization term which ensures the consistent dissipation of energy. 
Roughly speaking, the system (|l.lj) - (|1.3j) is governed by the coupling of the hydrodynamic fluid flow 
and the bending elastic properties of the vesicle membrane. The resulting membrane configuration 
and the flow field reflect the competition and the coupling of the kinetic energy and membrane 
elastic energies. 

For the system ()l.ip - (ll.3|) subjecting to no-slip boundary condition for the velocity field and 
Dirichlet boundary condition for the phase function, Du, Li and Liu in [1] obtained that there 
exists global weak solutions via the modified Galerkin argument, and there holds basic energy 
inequality 



|(i||„(.,t)||i a + 25(0(-, t))) +H|VuM)||i 2 + 7 ||^M|| 2 2 =0, V t > 0. (1.12) 



Moreover, the authors also proved that weak solution is unique under an additional regularity 
assumption u £ L 8 (0, T; L A (Q)). Recently, local in time existence and uniqueness of strong solution 
to the system (|l.ip - (|1.3p have been established in [Tj5], and under the assumptions that the initial 
data and the quantity (|fi| + a) 2 are sufficiently small, the authors proved existence of almost 
global strong solutions. Note that they have to restrict the working space with proper limited 
regularity due to some compatibility conditions at the boundary which is required in the fixed 
point strategy. Very recently, Wu and Xu [33] considered the system (|l.l[) - (|1.3j) with initial 
conditions (jl.4p and periodic boundary conditions (jl.5p to avoid troubles caused by the boundary 
terms when performing integration by parts. They proved that, for any given initial data (uq, </>o) £ 
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Hp er (Q) x Hp er (Q), there exists a positive time T such that the system (|l.ip - (|1.5p admits a unique 
smooth solution (u, <fr) satisfying 

u 6 C([0, T],H* er (Q)) n L 2 (0, T; i? 2 er (Q)) n H^O, T; L 2 per (Q)), ^ 
4>eC([0,T],H^ r (Q))nL\0,T;H^ r (Q))nH\0,T;H^ r (Q)). 

Moreover, if the viscosity /i is assumed to be properly large, then the system (|1 . If) — (|1 .5|) admits a 
unique global strong solution that is uniformly bounded in H^ er x H^ er on [0, oo). However, as for 
the well-known Navier-Stokes equations, an outstanding open problem is whether or not smooth 
solution of (|l.lj) - (|1.5[) on [0,T) will lead to a singularity at the time t = T. 

For the Navier-Stokes equations, some results were obtained in early by Prodi [35], Serrin [23] 
and Giga [12], they proved that if 

f T 3 2 

/ ||ti(-,f)||?, dt < oo with - + - = 1, 3<p<oo, (1.14) 
Jo P Q 

then the smooth solution u can be extended past the time T, while the limit case p = 3 was 

proved by Escauriaza et al. [10]. In 1995, Beirao da Veiga [2] established similar criterion for the 

derivative of the solution, i.e., (|1.14[) can be replaced by the following condition: 

f T 3 2 

J ||Vw(-,f)||^ dt < oo with - + - = 2, 3 <p < oo. (1.15) 

In 1984, Beale, Kato and Majda in their pioneer work [T] showed that if the smooth solution u 
blows up at the time t = T, then 

i-T 

IK-,t)||Loo dt = oo, (1.16) 



where u) = V x u is the vorticity of the velocity field. Later, Kozono and Taniuchi [16] and Konozo, 
Ogawa and Taniuchi [15] refined the criterion (|1.16p to 

\\u(;t)\\ B MO dt = oo and / \\uj{-.t)\\g dt = oo, (1.17) 



respectively, where BMO is the space of Bounded Mean Oscillation and B^ ^ is the homogeneous 
Besov spaces. Recently, Fan et al. [TT] and Guo and Gala [T3] improved the above criteria to the 
following two logarithmic type criteria: 



T IK-,*)IIbo 



o ^l + ln(e+|| U ;(-,t)|| i 30 o 



dt = oo (1.18) 



and 

\H;t)\ 



„ r dt = oo. (1.19) 

When the phase function <j> is considered, similar regularity criteria as (|1.14l) and (I1.15[) for the 
system ip — (|1 .3)) have been established in [25] . The first author of the present paper in [2S] 
obtained that the Beale-Kato-Majda criterion (|1.16[) still holds for the system (|l.ip ~ (|1.5p . 
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Due to the lack of global well-poedness theory of the system f|l.l[) (|1.5l) . the investigations of 
blow-up criteria of local smooth solution are very important ways to understand the properties of 
solutions. Motivated by the above results, the purpose of this paper is to establish the blow-up 
criteria for the system (ll.ip - (|1.5|) in term of the norm of the homogeneous Besov space. The main 
results of this paper are as follows: 

Theorem 1.1 For (uq, 4>q) G Hp er (Q) x H® er (Q) with V -uq — 0. Let T* be the maximal existence 
time such that the system (|l.ip - (ll.5l) has a unique smooth solution (u,(p) on [0,T*). If < oo, 
then 



IK* 



I B° 



-.dt = oo, (1.20) 



I B° 



.^l + ln( e +|M-,t) 
where uj = V x u is the vorticity field. In particular, 

limsup |K,£)||£o =oo. 

t/*T, ™'™ 

Theorem 1.2 For (uo,0o) G Hp er (Q) x HperiQ) with V-uq — 0. LetT* be the maximal existence 
time such that the system (II .![) - (II .51) has a unique smooth solution (u,<p) on [0,T*). If < oo, 
then 



IK,*)II| 



In particular, 



1 + ln(e+ 



limsup \\u}{-,t)\\g- 



-dt = oo. (1.21) 



Remark 1.1 Theorems 11.11 and 11.21 are still true, if we replace the vorticity ui by Vu in (|1.20[) 

and (|1.21|) . due to the boundedness of Riesz transforms in B^ p ^\q) and B^% er \Q)- For the 
definitions of these spaces, see Section 2. 

Remark 1.2 Since L°°(Q) B^ p ^\q), the result (TOO]) improves the Beale-Kato-Majda blow- 
up criterion in [25]. 

Remark 1.3 Observe that Vw G Bo^^ er \Q) is equivalent to u G B^ P ^(Q) and the Sobelev 
embedding L per (Q) B^^ er \Q). Therefore, Theorem 11.21 implies that if T* < oo, then 

,T. \H;t)\\% 
(<) J l + ln(e + \\u(-M^J dt = °' 

W Jo l+He+\\Vu(;t)\\L*) ' 

The rest of the paper is arranged as follows. In Section 2, we recall the Littlewood-Palcy 
decomposition and definition of the homogeneous Besov spaces, and review some crucial lemmas. 
In Section 3, we establish the bound of ||VA0|| L 2, which enables us to derive some specific higher- 
order energy estimates. In Section 4, we present the proof of Theorem ll.il Section 5 is devoted to 
the proof of Theorem 11.21 
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2 Preliminaries 



We first recall the Littlewood-Paley decomposition. Let iS(R 3 ) be the Schwartz class of rapidly 
decreasing function, and 6>'(]R 3 ) be its dual. Given / G S(IR 3 ), the Fourier transform of it, F(f) 
f . is defined by 

For any given g E <S(M 3 ), the inverse Fourier transform T~ x g = g is defined by 

F-Hg^x) = g(x) = / <?(£K*'« d£. 

Let Pi = {£ e M 3 , |C| < f } and £> 2 = {£ e R 3 , f < |£| < f }• Choose two non-negative radial 
functions </>, tp £ 5(R 3 ) supported, respectively, in 2?i and X>2 such that 

^(0 + $>(2-^) = l, £eK 3 , 

5>(2-'f) = l, fe« 3 \{0}. 

iez 

Let /i = J 7 ^ 1 ^ and /i = J 7-1 ^. Then we define the dyadic blocks Aj and Sj as follows: 

A,/ = 0(2-''D)/ = 2 3 ^ f h(Vy)f(x - y) dy, 

Sjf = ^(2- j D)f = 2 3 ^ f h(Vy)f(x - y) dy, 

where D = {D\, D2, -D3) and Dj — i^ 1 d Xj (i 2 = —1). The set {Aj, Sj}j^z is called the Littlewood- 
Paley decomposition. Formally, Aj = Sj —Sj-i is a frequency projection to the annulus {|£| ~ 2 J }, 
and Sj = Ysk<j-i ^ k i s a frequency projection to the ball {|£| < 2 J }. For more details, please 
refer to |18j . 

Next we recall the definition of homogeneous Besov spaces. Let 5^(R 3 ) be the space of tem- 
perate distributions / such that 

lim Sjf = in S'(M. 3 ). 

j^-00 

For sGM and (p, q) <E [1, 00] x [1, 00], the homogeneous Besov space Bp* g (K 3 ) is defined by 

B s p , q m = {/ G S' h (R 3 ) : Il/H^ < 00}, 

where 

= f(E ieZ 2^||A,-/||i P ) 1/9 for 1<?<oo, 
[sup jeZ 2J s ||Aj/|| L p for q = oo. 

It is well-known that if either s < | or s = | and g = 1, then (B^ ? (]R 3 ), || • ||^ s ) is a Banach space. 
In particular, when p — q — 2, we get the homogeneous Sobolev space H S (M. 3 ) = -B|.2(^ 3 ) which 
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is endowed the equivalent norm WfWjjs = ||(-A) s / 2 /|| L 2. The notation H S (R 3 ) is the standard 
inhomogeneous Sobolev spaces which is endowed the standard norm ||/||ir<> = || (— A) s / 2 /|| i 2 + 

\\f\\L>- 

We also need to introduce some well-established functional settings for periodic problems: For 
1 < r < oo, we denote by 

L r per (Q) ■= {u e L r (R 3 ) | u(x + et) = u(x)} 
endowed the usual norm || • \\l t - For an integer m > 0, we denote by 

H™ r {Q) := {u E H m {R 3 ) I u(x + (h) = u{x)} 
endowed with the usual norm For s 6 R and (p, q) £ [1, oo] x [1, oo], we denote by 

B'^HQ) = { u e B S P JM. 3 ) : I u(x + a) = u(x)} 
associated with the norm || ■ || f, B . 

Before ending this section, we state some well-known inequalities. The first one comes from 

Lemma 2.1 ( 14 ) For s > 1, we have 

l|V s (/ 5 ) - fV s g\\ LP < C(\\Vf\\ LPl W^-'gW^ + IIV/H^llflllw) (2-1) 
with 1 < w, oi , no < oo such that - = — + — = — + —. 

- fl IX! ft ~* p pi qi P2 <J2 

The second one can be found in [15] and the proof follows from the Littlewood-Paley decom- 
position. 

Lemma 2.2 ([15]) For all f E iJ^^M 3 ) with s > §, we have 

ll/lk- < C(l + Il/H^ ln^ 2 (e + II/Hh-0)- (2-2) 
The last one comes from [3D], see also [T3"] . 
Lemma 2.3 ([20] [15]) For all f E i7 1 (M 3 ) n ^"^(K 3 ), we have 

ll/IU*<C||/||V| v ||/||^ 2 . (2.3) 

3 The bound of ||VA^|| £ 2 

By the basic energy estimate (|1.12p . we can easily get the following uniform estimates (cf. [4l [25]): 

\\u(-,t)\\ L 2 + \\(f>(-,t)\\ H 2 < C for all t>0, (3.1) 

l (v\\Vu(-,t)\\h+l\\^(-,t)\\h)dt<C, (3.2) 

where C is a constant depending only on ||mo||l 2 j ll'/'olli? 2 an d coefficients of the system except the 
viscosity /i. 
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Lemma 3.1 Assume that (uq, <po) 6 Hp er (Q) x H® er {Q) with V • uq = 0. For any smooth solution 
(u,4>) to the system (|l.l[) - (jl.5l) . we have 

sup ||VA0(-,t)|| L2 <C (3.3) 

0<t<T 

for any < T < oo, where C is a constant depending only on \\uq\\iji, ||0o||_ff 3 > T and coefficients 
of the system. 

Proof. Taking A on (|1.3[) . multiplying the resultant by — A 2 </>, and integrating over Q, we obtain 
--||VA^||i 2 = - / V • (u • V<j>) • VA 2 ^ - 7 / V x --VA 2 <^:=/ 1 +/ 2 . (3.4) 

For Ji, by using the interpolation inequality ||V 2 (/>|j 2 3 < C|| V 2 </>||£ 2 1| VA0||x,2, we can infer from 
(|3~T]) that 

/i < ^||VA 2 0||i 2 + C||V« • W||| 2 +C\\u- V 2 0||i 2 

< ^||VA 2 ^|| 2 2 + C(||V U ||i 2 ||V0|| 2 . + h||ie||V 2 0||i3) 

< ^£||VA 2 0||| 2 + C||V U || 2 2 (||VA0|| 2 2 + 1) + C||V u || 2 2 ||V 2 0|| L2 ||VA^|| i2 

< ^£|| V A 2 0||i 2 + C||V U ||| 2 (||VA^||| 2 + 1). (3.5) 
For J2, since A(4>) and B{<p) are functions depending only on time, by we obtain 



7 / V 

Q 



kg{4>) + Mi(A(^) - a) + M 2 (B(0) - /3)/W>) 



■ VA 2 <Mr 



[ VA/(0) • VA 2 0da; - -J / V[(3</> 2 - l)/(<£)] • VA 2 0dx 



fc 7 / VA/(0) 

'Q 

M 2l (B(cf>) -(3) f V/(0) • VA 2 </)dx 



:-/ 2 i+/ 2 2 +/23- (3.6) 
Note that f((j>) = —eA<fi + -j(<^ 2 — l)</>, by (13.11) . we can estimate In (i — 1, 2, 3) as follows: 

hi = -fc£7||VA 2 </>|| 2 2 + ^ f VA(0 3 - (f) ■ VA 2 ct)dx 
- - — i||VA 2 0|| 2 2 + C||VA(0 3 -0)|| 2 2 



7fc£7 n vv A 2 j 1 1 2 , rif\\A\\4 IIV7AJ.II2 1 ||i||2 ||V7^I|2 II A ^l|2 



< — g-^nvA 2 0ii| 2 + c(u\\u\^^\\h + 

+ \\V<t>\\ie + ||VA0|| 

L 2 ) 

< _Z^I||VA 2 0|| 2 2 + C(||VA0||| 2 + 1); (3.7) 



/ 22 = / ■ VA 2 <f>dx ~m (3^ 2 - l)V/(0) • VA^cfe 

e Jo £ Jo 



< _£I|| VA 2 ||2 2 + W/(0)||| 2 + ||(30 2 - 1)V/(0)|| 



< 



||VA 2 0|| 2 2 + C(|| 0V0A0||i» + ||^ 2 (0 2 - 1)V^|| 
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+ ||(30 2 - 1)VA0|| 2 2 + ||(3^ 2 - 1)V(^ 3 - 0)|| 2 2 ) 

< ^Z||VA 2 0|| 2 2 + c(||0|| 2 oo ||V0|| 2 e||A0|| 2 3 + II^Hl^H^ 2 - l||i^||V^||| 2 
+ ||3^ 2 - 1|| 2 .||VA^|| 2 2 + ||30 2 - 1||U|V^|| 2 2 ) 

< ^||VA 2 0||| 2 + C(||VA0|| 2 2 + 1); (3.8) 

Jaa = M 2 £f(B(<f>) -13) f VA0 • VA 2 0dx - M ^ B ^ " ® f V (0 3 - 0) • VA 2 (f>dx 
Jq e Jq 

< ^2|| V A 2 0|| 2 2 + C(B(d>) - /3) 2 (||VA^|| 2 2 + ||V(0 3 - ^)|| 2 2 ) 

< ^2|| V A 2 0|| 2 2 + c(||V^||i 2 + ||0 2 - + l) (||VA0|| 2 2 + ||0 2 - 1|| 2 » ||V0||| 2 ) 

< ^2|| V A 2 0|| 2 2 + C(||VA0|| 2 2 + 1). (3.9) 



From (|3I71)-(1SU), we get 

/ 2 <-^||VA 2 ^|| 2 2 +C(||VA0 i.. :i). (8.KI) 

Combining the above estimates (|3.5[) and (13.101) . we obtain 



|||VA0||| 2 + fe 7 ||VA 2 0||| 2 < C(||V U || 2 2 + 1)(||VA0||| 2 + 1). (3.11) 
The Gronwall's equality yields that 

||VA0(t)||| 2 < ||VA0o||i»exp(cjr (||V«(r)||| 3 + l)dr). (3.12) 

The estimate (|3. 121) with (|3.2j) imply (|3.3[) immediately. We complete the proof of Lemma [3. II □ 
By (|3Tj) and (|3T3j) . for any < T < oo, we obtain 

sup U(;t)\\ H 3 <C. (3.13) 
0<t<T 

By the Sobolev embedding H 2 er (Q) ^ L™ r {Q), (|37T3]) yields that 

sup ||V0|| L = < C*. (3.14) 

0<t<T 



This result will be used frequently in the proofs of Theorems 11.11 and 11.2 



4 The proof of Theorem 11.11 



We argue Theorem 11.11 by contradiction. Assume that the result (|1 .20|) is not true, which means 
that there exists a constant M > such that 



r 

Jo 



IM-,*)IIb° 



-.dt < M. 



Under the condition (|4.1[) . if we can prove that 

IimSU P (||u(-,t)|| ff 3 + \\<P(;t)\\ H e) < C 

tyr, 



(4.1) 



(4.2) 



holds for some constant C depending only onuo, 00: M, T* and coefficients of the system (|1.1[| - 
(|1.5I) . then we can extend the solution (u,<p) beyond the time t = T„, which leads to the contra- 
diction. Therefore, it suffices to show that under the condition (I4.1[) . we get 



Taking the curl on ()1.1[) . we obtain 



SE. 



d t uj - (aAuj + u ■ Vw = w ■ Vu + V x (-ttV</>) 



(4.3) 



Multiplying ()4.3j) by uj and integrating over Q, we have 

' "-MH, + HIV.. 7, = 

IQ 



2 dt 



f f SE 

/tillVwII n = / w-Vu-iodx— I ——S7(f)-\7xu)dx, 
Jo Jo 6( P 



(4.4) 



where we have used the fact Jq u ■ Vw • uidx — due to V • u = 0. Since the Riesz operators are 
bounded in L 2 and Vu = (— A) _1 V(V x ui), we have ||Vu||ia < C||w||l2. This implies that 



w ■ Vw • udx 



<C||w||i-||V«|U a ||w|U»<C||w|| L =.||w|| 



Applying Young's inequality and (|3.14l) . we have 



SE 

— — Vq> ■ V x ujdx 



< J||V^||i 2 +C||^V0||| 2 

< J||V W ||i 2 +C||^|| 2 L2 ||V0||| o 
<fl|Vu;||£ 2 + C|fil 2 . 



Taking (|4"3)l and into ([O]) . we obtain 



+ f l|v w ||i 2 <c(|M| L » + i)(|Hii 2 + ||^' 



(4.5) 



(4.6) 



(4.7) 



On the other hand, after integration by parts, we obtain from (|1.11[> that 
1 d_ SE f ^S_E SE 



10 



^ [kg(4>) + Mx{A{<t>) -a)+ M 2 (B(4>) - /?)/(</>)] • ^dx 



+ M^M) I f{4>) ■ S -^dx + M 2 (B(<P) -13) [ • S -^dx 



' dt Jq S(f) ' J Q dt 

■■= Ji + J2 + J3 + Ji + 'h- (4.8) 

Noticing from (| 1 . 3[) that 

II^IU»<c(||ti-v0|U 3 + ||^|| L! i) 

,SE, 

,||V0||£»+| 

,6E, 



<c(H| ia ||V0|U- + 11^11^ 



||V^|| L2 < C(||V« • V0|| i2 + ||u • V 2 0|| i2 + ||vS La 

< C(||V«||i»||V0||i- + || W || L 3||V 2 0|| L 6 + ||V— \\ L 2 



::f(||A— \\ L2 + \\Vu\\ L2 + \\ — 



dt \\»<c(\\A^\\ L , + \\^-mv) 

< c(\\Au ■ V0|| L2 + || V« ■ V 2 0|| i2 + ||u • VA0|| L2 + ||A^|| i2 + ||^(0 3 - <j>)\\ L z) 

< c(\\V<j>\\ L ~\\Au\\ L , + ||V U || L3 ||V 2 0|| L6 + hmiva^ii^ + ||A^|U 2 + ||^|U= 



SE „ „<5.E 



<c(||A U || i2 + ||A— + + I 



Then we can estimate (i = 1, 2, 3, 4, 5) as follows: For J\, we can further split it into the following 
two terms: 

J 1= ke f A^ • A^gdx - * / |A(0 3 - 0) • S -£-dx := J u + J 12 . (4.9) 
Jq at 09 e J Q dt S<p 



By using Leibniz's rule, (|1.3|) yields that 

5E f SE 

J n = -kej\\A— 1|| 2 -ke A(u ■ V4>) ■ A—dx 

d(p Jq d(f) 

<-^l|A^||! 2 +C||A( w .V0)|| 2 2 
9fce7 „ A SE 



||A— 1| 2 2 + a(||A« • V0||| 2 + 2||V« • V 2 </>|| 2 2 + \\u ■ VA0||| 2 ) 



10 

11 



< 



< 



9fce7 ,, . SE ,,, 
'||A- "- 



10 

9kej 
10 



i .^c(||V^||l oo ||A M ||| 2 + ||V u ||! 3 ||V 2 0|| 
C||A U ||| 2 + C(||V U || 2 L2 + 1), 



2 



||«||i,||VA^|| 



2 

L 2 



(4.10) 



k f d ^3 



■ A—dx 



<§|| A «|| i , + c( M t.|«li, + ||| 

<§iiAa J . + c(n& I . + i 



Hence, we infer from (14.10)) and (|4.11[) that 

Ji<-^||A^||| 2 +a||A M ||i 2 + c(||V u || 
Similarly, we can estimate J2, J3, J4 and J5 as follows: 



dm SE 



6fc /" , A , 90 SE , 6fc 



2 

L 2 



90 <5£ fc /• 2 



(4.11) 



(4.12) 



-1) 



9/(0) 



< 



C(||0||^||A0|| L6 ||^|| L3 + ||0||| TC ||0 2 -1|| 
< C 



£-11-57 IIL 2 + II0 2 -1|U' 



0t 



1 



-rrdx 



£= 



|L 2 



^ llL2 + l|V W lli2 



W) N,^., 



< c(||Au|| ia + l|A^|| i2 + ||Vtt|| La + H^IU 2 + l) II^IIl 2 



50 



< ^W^Wh + CII** + ^(llVn||i 3 + + 1 



(4.13) 



J3 <C 



<L4(0) 




< c 


f 




dt 






i Q Tt dx 





SE 



lk.<c(||Sii, + i) 



L 2 / ' 



(4.14) 



< c(||V0iu 

<C(||A^ 
10 11 i 



ll/WOIMI^IU 2 
L2 ||V^|| i2 + ||0 3 - 011^11^11^) ||/(0) \\ L2 1| S ^\\ L2 

SE n »„ 11 11^,1 ^n^n 
"£ 2 + l|Vw|| i2 + || — || i2 + lj|| — || i2 

<^||A^||i 1 +(7(|[V* + ||S|i, + l) > 



(4.15) 



12 



j*<cm)-m^\\»\\ s 4\\* 



11 at 

<c(||V0||| 2 + ||0 2 -l||i 2 )|| 



dfW SE 



< c(\\Au\\ L * + l|A^|| L2 + ||Vu|| La + II^IU. + l) ||^|| £ > 

k£ luA SE u2 , ^MA„.||2 , „/ M __ ||2 ,,1^1,2 



- To i|a ^ IIl2 + c ll Au ^ 2 + c {\^\\h + II ^ Ilk + !j • ( 4 - 16 ) 

Taking (|4T2|) - (|4T6|) into flUJ) , by using the fact that ||Aw|| L 2 < C||Vw|| L 2, we conclude that 

+ fe 7l|A^||| 2 < <5||VHIi> + c(|Mli, + ||^||| 2 + l), (4.17) 

where C is a constant depending only on ||tto||jyi, ||0o||_f/ 3 j and coefficients of the system due to 
the estimate (12.31). 



Set 

= JL 

V 2(7 

Then multiplying (|4. 1 T[) by 77, adding (14.71) together, applying Lemma \2. 2 1 with s — 3, we obtain 
d /„ „, \ „„ „, , „ A SE, 



2 
L 2 



^(ikll^ + ^l^lll^+MllVHI^ + ^IIA^II 

<^(lHU- + i)(|Hli a +T 7 ||S|i, + i) 



< c(l + Hl^^l + lnCe+HI^)) (|M& + vW^Wh + 1 

< C(\\w\\l 2+V \&\l 2 + l)+C - l|W|lj? 7°°„ Jl + ln(e+HlBO 

x Vla(e+H|H»)(|klli a +»7||^lli a + 1 

<ff(llHli» + "11^11^ + 1 

+ g , = 7'°°,, ln(e+||V 3 ui| L2 )(|| W H| 2+?? |l — 1|| 2 + 1), (4.18) 

./l + ln e + H jjo ) V ^ 7 

y oo,oo 

where C is a constant which may depend on r\. 

By the condition (jl.20j) . one concludes that for any small constant a > 0, there exists T < T 
such that 

IMIso 

-.dt < a. (4.19) 



2 



To /l + ln(e+|H|BO ) 

For any T < t < T, we set 

**(«):= sup (||VA^t)||| 2 +^||A^(t)||£ 2 Y (4.20) 

T <T<t V 0© / 
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where fj is a determined constant which specified later. Applying Gronwall's inequality to (|4.18l) 
in the time interval [7b, i], one has 

\Mt)\\h + v\\^(t)\\h 



<C exp( Cds + C\n(e + h(t)) dr 
K Jt Jt /l + ln(e+ INIIbo,^) 



< C exp (C(t - T ) + Co ln(e + ft(t))) 



< C {e + h(t)) 2C ° , (4.21) 

where Co = ||u>(7b)||| 2 + 7 7l|^(7b)|| 2 2 is a positive constant depending on T . 

Now we are in a position to derive higher order energy estimates of the solution. Taking VA 
on (|1.1|) . multiplying VAit and integrating over Q, we obtain 



2~dl 



f f SE 

||VAu||? 2 +^||A 2 u||? 2 = - / VA(h • Vu) ■ VAudx + / VA( — V0) ■ VAudx 

' ' " Jo Jo S( P 



>Q J Q 
h+h- (4.22) 



Since V • U = 0, I\ can be rewritten as 

h = - [VA(u • Vu) - u ■ VVAu] • VAudx. 
Jq 

By using Lemma 12.11 we can estimate I± as follows: 

h < C||VA(u • Vu) - u ■ VVAu|| L 4/3||VAu|| L 4 < C|| Vu|| L 2 1| VAu||| 4 
< C||V M ||^ 6 ||A 2 U ||^ 6 < |||A 2 M || 2 2 + C\\Vu\\% 

<%\\A 2 u\\i 2 +C\\w\\l\, (4.23) 
where we have used the Gagliardo-Nirenberg inequality: 

||VAu|U« <C||Vu||^ 12 ||A 2 u ||^ 2 /12 . 
For I2, after integration by parts, by using (|3.3p and p. 141) . we obtain 

h = - I A{^-V<P)A 2 udx 
Jq o<P 

<|||A 2 U || 2 2 + C||A(^V0)|| 2 2 

< |||A 2 u|| 2 2 + c(||A^|| 2 2 + 2||V^V 2 ^|| 2 2 + II^VA0|| 2 2 ) 



< -IIA 



2„,l|2 



U 



L 2 



C(||A^|| 2 2 ||V0|| 2 , + ||V^|| 2 3 ||V 2 0|| 2 6 + ||^|| 2 .||VA0|| 2 2 



< ^\\A 2 u\\h + cf||A^||| 2 + ||^|| 2 L2 + l). (4.24) 
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Combining (|4T2l2l - P~2"4l . we deduce that 
d . Il9 3/i 



|||VAu||| 2 + f ||A 2 W || 2 2 < + HA^Hi, + f-^Wh + 1 



SE, 



(4.25) 



To obtain the desired estimates for (f>, we start from that 
1 d „ . SE , 



2 dt 



||A- 



2 
L- 



d . SE ^SE , 



= y ^A[Atf(0 + Mi(A^) - a) + M 2 (B(0) - (3)f(<f>)] ■ A^dx 



-Q^_A[kg((f>) + M 2 (B((f>) - /3)f ((/>)] ■ A—dx 



: J §1 [%W + M 2 (B(<f>) ~ ■ A 2 ^dx 

, f 9 . ., ,, . n SE k 
-kj Q -A m .^- dx + 7i 

+ M 2 jB{4>) J f{<p) ■ A 2 ^-dx + M 2 (B((f>) - ft) J ^-f^.A^dx 



:= J1 + J2 + J3 + Ja- 



(4.26) 



Let us estimate the terms Ji (i — 1, 2, 3, 4) one by one. For Ji, we divide it into the following two 
parts: 



11 "t" Jl2 



(4.27) 



For Jn, by using Leibniz's rule, we deduce from (| 1 .3[) that 

Ju = -fce 7 ||A 2 ^||| 2 - fee / A 2 (u • V0) • A 2 ^-dx 



< 



< 



^\\A 2 ^\\l 2 + C\\A 2 (u-^)\\h 

A 2 ^||| 2 + c(\\A 2 u ■ V0|| 2 2 + 4|| VAu ■ V 2 0|| 2 2 

+ 6|| Au ■ VA<£||| 2 + 4||Vu • S7 2 A(t)\\ 2 L2 + \\u ■ VA 2 0||| 2 ) 

15fc£7 
16 11 8<t> 

+ ||A U || 2 6 ||VA0|| 2 3 + ||V U || 2 TC ||A 2 0|| 2 2 + \\u\\ 2 L3 \\V 5 cp\\ 2 La 



— A- > ^|| 2 2 +C(||V0|| 2 =O ||A 2 W || 2 2 + ||VA W || 2 3 ||V 2 



\L' : > 



< 



< 



15fc£ 7 |, A 2^||2 



16 



\\i 2 +C\\A 2 u\\i 2 +C 
5E, 



1)(||VA W || 2 2 + 1) 



(||Vz 1 || 2 2 + l)(||A— 1| 2 2 + 1) 



16 

•C(||V«|| 



2 

L- 



L 2 -tC\\A 2 u 
6E 



2 , ||2 
L 2 



\h 



1)(||VA M || 2 2 



A- 1,2 



(4.28) 
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where we have used the facts || A 2 0||2 2 < C(|| || 2 2 +l) and || V 5 0|| 2 6 < C\\V 6 <j>\\ 2 L2 < C(||Af|||| 2 + 
1). For J12, it clear that 

f k f 5 a/,3 ^ a2^ j 6fc [ d(\V(f>\ 2 <f>) a2 5E , 



3k r d(<j) 2 A(j)) A 2$E_ dx+ k f A 2^_ dx 
e Jq dt 84> e J Q dt 5<j) 

12k f ^^dej) a2 5E , 6fc [ .„ ll2 d<t> a2 5E , 
e J Q dt 5cj> e J Q dt 8<j> 

St- ■ A 2 —dx / <i 2 A-r- • A —-dx 



e Jq dt 8<j) £ Jq dt 



+ -/Af .A'gd* (4.29) 



Similarly, we can estimate the terms Ji2i (« = 1, 2, 3, 4, 5) as follows: 

^^l|A 2 fll| 2 ^IIWvf||| 2 

<^||A 2 ^||| 2 + C||^||V^||V^||| 2 



< 5 l|A2 f lli2 + C ( l|A f »' a + l|Vu|li2 + "f + 0' (4 ' 30) 

/ ^ — ||A 2 f || 2 ,-q|Vo||Mit -f. 



16 11 ty Mi V" 50 



Ji23<^||A 2 ^||| 2 + C||^||A^|| 2 6 ||^||| 3 

^li^lii- + ^(n«-v* + ii^lii.) 

<^l|A^||i, + c(|N|i,|^||i. + ||Ag|,i, + ||g„i,) 

< ^IIA 2 ^ IIL + C(||V W ||| 2 + ||A^||| 2 + ||^|| 2 2 ), (4.32) 



Ji 24 <^||A 2 ^|| 2 2 + q|0|| 2 »||A^|| 2 2 

< ^W^Wh + c( l|A(u • *4>)\\h + l|A^||| 2 



50 

< ^H A2 ^lli^ + c(l|Vtf||i~||Au||£ a + ||V U || 2 3 ||V 2 0|| 2 6 



+ |M| 2 .||VA0|| 2 2 + ||A— 1|| 2 ) 
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< ^||A 2 S|| 2 + c(||A«||i, + HAglli, + l) , (4.33) 



16 11 8(t>^ V" 11 8<t> 



f < k£1 \\A 26E \\ 2 4-r\\A d<j} \\ 2 



< ^W^Wb + c(||Au||i, + ||A^||£ a + l). (4.34) 
Putting estimates (|4.30p - (|4.34|l together, we obtain from (|4.29|) that 

Ju < ^p||A 2 ^||| 2 + c(||Au||| 2 + ||A^||| 2 + ll^Hi, + l). (4.35) 
Since ||Au||| 2 < C(||VAu|| 2 2 + 1), we obtain from (|4~551) . (I4T2T)) and (|4~25| that 
J^-^IIA^IIi. + GHA^Hi, 

+ c(\\Vu\\l 2 + ||^|||. + l) (||VA«||£ a + ||A^[[£ a + l) . (4.36) 
The estimates of Ji (i = 2, 3, 4) can be proceeded as that of Ji (i = 2, 3, 4, 5), thus we have 



J 2 = §/ W )g.A^4/(3^-l) 



dt 5$ 

6k f , . , dd> . o SE . 6k f , 2 . , 2 , . dip A? SE , 
(pA<j)— - ■ A —dx + <t> {4> - l)^f ■ A —dx 



e Jq dt S(f> s 3 Jq dt 



3 / ^ 2 - D 

£ JQ 



at <w 



< 



:| ^lli* + c(MU\Anl e \\^\\ 2 L3 + H\\U\4> 2 ill 2 



+ m 2 -nu\^\\h) 



^Y" A ^^ + c UI^II- + ll A W"--"at"- 



< Y lA^IlL + c(||VAu||£, + IIA^IH. + Q|| 2 + l), (4.37) 



i3<C|^«|||/ W || L2 ||A 2 ^|U 2 



< YHA 2 ^||i 2 + c(j|v</,Mlv^|| i2 + w - ^ii^ii^j wmwi* 

< ^HA 2 f \\h + cfllV* + ||A^||| 2 + ||f ||| 2 + lV (4.38) 



j 4 <a TO )-^|||^M|u 2 ||A 2 ^|U 2 



< 



^l|A^||! 2 + c(||v^||! 2 + ||0 2 -i||i 2 ) 2 ||^[ 



\h 
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<^l|A 2 Sl!l 2 + ^l!A#ll! 2 + ll#ll! 2 



< 



at"" "at 



' l|A u lli2 + C (" VA * + » A D i2 + Ql* 2 + ^ 



Taking (jOS|) - (|Qg)| into (|4T26]) . we conclude that 

|<B. + ^|A^li.<<S||A'.f 1 , 

+ C(||V«|||, + |^|||, + l)(||VA»||i, + IIA^Hi, + l), (4.40) 

where C is a constant depending only on HitoH/fi, ||0o||if 3 j T* and coefficients of the system due to 
the estimate (|3.3[) . 
Set 

- = JL 
V 2C' 

Multiplying (|4.40|) by fj and adding the resultant to (|4.25f) . we obtain 



(||VA«||i, + A^||| 2 ) + „|| A 2 M || 2 2 + A 2 ^||| 2 



< c(||Vu||l 2 + ||^||| 2 + l) (|Mlli + l|VAn||l 2 + flW^Wh + e) • (4.41) 



It follows from (|4~2Tj) that 

j t {e + h{t)) < c(||Vu||| 2 + H^lli, + l) ((e + h{t)f^ + h(t) + e) . (4.42) 



Choosing a small enough such that lACa < 1, we get 



j t (e + h(t)) < c(||Vu|| 2 2 + ||^|| 2 2 + l) (e + h(t)). (4.43) 



Applying Gronwall's inequality leads to 

h(t) < (e + h(T ))exp (c J (\\Vuf L2 + ||^|| 2 L2 + l)dr 



(4.44) 



This combines with the basic energy (|1.12p yield the boundness of h(t) on the time interval [Tb,T] 
Since it is easy to verify that 



r pi 

\\<f>\\%e <C(||A— || 2 2 + 1), 



we finally obtain from (|4.41j) that 



SUP (\\u(;t)\\ H 3 + U(;t)\\ H e) <C. 
T <t<T 



This completes the proof of Theorem 11.11 
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5 The proof of Theorem 11.21 



Similarly we prove Theorem 11.21 by contradiction. It suffices to prove that if 

\H;t)\\U 

dt<M< oo, 



o 1 + ln(c + ||w(-,*)IU-x_) 



(5.1) 



then 



limSUp(]K-,t)|| H 3 + \\(P{;t)\\ H e) < C 



(5.2) 



for some constant depending only on uq, 4>q, M, T* and coefficients of the system (|l.ip - (|1.5l) . 
Multiplying (|4.3p by ui and integrating over Q, we have 



1 d f f SE 

^tIMI? 2 + mIIVwII ? 2 = / w ■ Vu • uidx - / — Vr/> • V x ujdx, 

2 dt J Q J Q 5<p 

Since ||Vtt[[ia < C||a;||x,2, by using Lemma 12.31 we obtain 

w ■ Vu ■ ujdx < C||w||| 4 ||Vu|| L 2 < C|M|| 4 |M| L 2 
<C|HIb-%„IHIl»I|VHIl» 

^fllvHIia + cHII- Hi!- 

The second term on the right-hand side of (|5.3[) can be estimated the same as (|4.6I) : 



(5.3) 



SE. 



-V</> ■ V x wdx 



Taking ([5TI|) and ([53)1 into ([Of . we obtain 
5lHli- + ^HV«||i-<C7 



<Sllv»[ii- + cf||ftii,. 



i)(Hli- + llf 



(5.4) 



(5.5) 



(5.6) 



The estimate for jt can be proceeded the same as that in the proof of Theorem 11.11 thus we also 
get 1|4.17[) . Multiplying (|4.17[) by r\ and adding (|5.6[) together, we obtain 



!(Nli* + »?ll^ll£») +Mlivn>|li» +ks^\\£^f L , 

SE 



,SE, 



SE, 



<c(||c|||_ loc + l)(|h||| 2 + ?7 ||^|| 2 L2 + l 



C 



< c 



1 + He+\Hs-]J 
l + ln(e+||w|| A -i ) 



hE 

Mh+VWT-Wh + I) l + ln(e+||VAu|| iB ) 



H\ 2 L 2 



SE 



L- 



l + ln(e+ ||VAu|| ia ) , (5.7) 
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where C is a constant which may depend on r/. 

By the condition (ll.2ip . one concludes that for any small constant a > 0, there exists To < T 
such that 

.T H">l| 2 5 +1 

dt < a. (5.8) 



b l+Me+||w||B-i ) 



For any T < t < T, we set 

up (||\^/(r! 7, f //|p' >£ 



M*):= sup (||VAu(r)||i 2+7 }||A^(r)||i 2 ), (5.9) 

T <r<t V 0(j> J 



where r) is a determined constant which specified later. Applying Gronwall's inequality to (|5.7p in 
the time interval [Tb,i], one has 

\\Vu(t)\\l 2+V \\ — (t)\\l 2 



t r nil- + 1 \ 

< Coexp(C(l + ln(e + fc(i))) / . „ -dr) 



< C exp (Ca(l + ln(e + h(t)))) 

<C (e + h(t)) 2Ca , (5.10) 

where Co = ||Vu(T )||| 2 + 7 ?ll7jf'(T'o)|||2 is a positive constant depending on T . 

The derivations of higher derivative estimates are analogously the proof of Theorem thus 
we safely omit it. This completes the proof of Theorem 11.21 
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